Frequency-dependent ratchet effect in superconducting films with a tilted washboard 

pinning potential 
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The influence of an ac current of arbitrary amplitude and frequency on the mixed-state dc-voltage- 
ac-drive ratchet response of a superconducting film with a dc current-tilted uniaxial cosine pinning 
potential at finite temperature is theoretically investigated. The results are obtained in the single- 
vortex approximation, i.e., for non-interacting vortices, within the frame of an exact solution of the 
appropriate Langevin equation in terms of a matrix continued fraction. Formulas for the dc voltage 
ratchet response and absorbed power in ac response are discussed as functions of ac current ampli- 
tude and frequency as well as dc current induced tilt in a wide range of corresponding dimensionless 
parameters. Special attention is paid to the physical interpretation of the obtained results in adi- 
abatic and high-frequency ratchet responses taking into account both running and localized states 
of the (ac+dc)-driven vortex motion in a washboard pinning potential. Our theoretical results are 
discussed in comparison with recent experimental work on the high-frequency ratchet response in 
nanostructured superconducting films [B. B. Jin et ai, Phys. Rev. B 81, 174505 (2010)]. 

PACS numbers: 74.25.F-, 74.25. Wx, 74.25. Qt, 74.40.De 



I. INTRODUCTION 

Within the last decade vortex ratchets, which exploit 
asymmetric vortex dynamics, have been attracting con- 
siderable attention 1 - 13 . In essence, the vortex ratchet is 
a system where the vortex can acquire a net motion in 
an asymmetric periodic pinning potential (PPP) in the 
presence of deterministic or stochastic forces with time 
averages of zero (for comprehensive reviews, seei-^). The 
asymmetry in the PPP refers to the current direction re- 
versal. There are essentially two different ways to realize 
such a pinning potential asymmetry. First, the spatial 
inversion symmetry of the PPP itself can be broken in- 
trinsically and involves some kind of periodic and asym- 
metric pinning potential, also known as a rocking ratchet. 
A second option is that an initially symmetric PPP, if 
externally biased, i.e., subjected to an additive constant 
driving force, results in an effective asymmetric pinning 
potential. This is called a tilted-potential ratchet or a 
tilting ratchet. Irrespective of the way to bring the asym- 
metry into a system, from a practical viewpoint, the com- 
mon feature of superconducting ratchets is their rectify- 
ing property: the application of an ac current to a super- 
conductor with an asymmetric PPP landscape can pro- 
duce vortex motion whose direction is determined only 
by the asymmetry of the pinning potential. 

A considerable amount of theoretical work about the 
general properties of different types of ratchet systems 
existsi^. Such ratchet systems range from the use 
of Josephson junctions in superconducting quantum in- 
terference devices (SQUIDs) and arrays^, to the use of 
one- or two-dimensional potential-energy ratchets allow- 
ing one to construct fiuxon pumps and lenses^, and drive 
fluxons out of superconducting samples^. Whereas the 
majority of ratchet proposals rely on single particles in- 



teracting with an external potential to produce the dc 
response, collective interactions between particles needed 
to produce dc transport have also been considered previ- 
ously 8 . At the same time in experiments, initially asym- 
metric PPPs have been used rather than dc-biased ones 
so far. For instance, the vortex lattice ratchet effect 
has been investigated in Nb films sputtered on arrays 
of nanometric Ni triangles, which produce the asymmet- 
ric PPP 9 . Similar effects were also discussed for YBCO 
films with antidotal. The voltage rectification in super- 
conducting Al films patterned with either asymmetric or 
symmetric antidots and subjected to an ac driving cur- 
rent, dc-biased in the symmetric case, has been experi- 
mentally observed in RefJ^. Among other experimental 
works on the vortex ratchet effect in nano-patterned su- 
perconductors, two substantial recent paper s 1213 should 
be mentioned. D. Perez de Lara et al*£ have investigated 
ratchet effects in thin Nb films grown on top of arrays of 
Ni nanotriangles subjected to an ac current with a fre- 
quency up to 10 kHz, so that effects observed in that 
worki 2 ^ were adiabatic, i.e., independent on ac frequency. 
Only recently, B. B. Jin et al. 13 have experimentally in- 
vestigated a very important issue in the vortex ratchet 
study, namely the frequency dependence of the dc voltage 
at large amplitudes of the ac driving force in a frequency 
range between 0.5 MHz and 2 GHz. As it was pointed 
out in RefJ^, ac frequencies were always lower than 1 
MHz in vortex ratchet measurements up to that work 13 
(see Refs. [14-20] therein). 

So far, a full temperature-dependent theoretical de- 
scription of the superconducting devices proposed in 
Refs. [9-13] is not available due to the complexity of the 
two-dimensional PPP used in these references. In par- 
ticular, the theoretical explanation of the experimental 
results of the vortex flow along the vortex channeling 
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directions in the above-mentioned structures is a diffi- 
cult problem. Below we propose to experimentally study 
ratchet properties on nanostructured thin-film supercon- 
ductor s 14 ' 15 with an uniaxial, i.e., washboard pinning 
potential (WPP). Uniaxiality of the proposed potential 
does not mean that the physics of vortex motion in a 
symmetric WPP, tilted by an external dc bias, becomes 
one-dimensional as the angular component of the moving 
force could simply be regarded as changing the strength 
of the WPP. Besides the two limiting cases of transver- 
sal and longitudinal geometry, when the vortices move 
across or along the WPP channels, respectively, one has 
to consider for all intermediate angles not only longitudi- 
nal, but also transverse ratchet effect^^, where the last 
appears due to the guided vortex motion along the WPP 
channels. 

In the present work, however, we will use the trans- 
verse geometry only to provide the reader with most intu- 
itive figure data and to simplify the subsequent analysis. 
Whereas the model allows one to obtain exact expressions 
for the magneto-resistivities at any intermediate current 
angles with regard to the WPP channels, general for- 
mulas for the responses have been provided recently^. 
In that worki£, a full and exact theoretical description 
of the nonlinear in current and temperature vortex dy- 
namics in the ratchet devices have been performed in the 
single- vortex approximation, i.e., for non- interacting vor- 
tices, within the framework of the Langevin equation. It 
should be noted that theoretical predictions of Ref.— lack 
experimental scrutiny so far, especially at microwave and 
GHz frequencies. For this reason, throughout this paper 
where applicable, we will be referring to the results of 
B. B. Jin et al^ as a most closely related experimen- 
tal work to our new tilted-ratchet results, though that 
worki 3 - is dealing with a rocking ratchet. 

In more detail, the mixed-state resistive response of su- 
perconducting films has been theoretically investigated in 
Ref.— in the high frequency and strong amplitude regime 
of the ac vortex transport in the presence of a dc bias, 
which invokes a definite tilt of the cosine pinning poten- 
tial, taking nonzero temperature fluctuations also into 
account. The exact solution of this nonlinear and time- 
dependent problem, obtained in Ref J£ in terms of a ma- 
trix continued fraction, included only general formulas 
for the ac and dc magnetoresistive responses. These will 
be elaborated in detail in this paper and applied to the 
study of both the dc ratchet electric field response and 
the absorbed ac power dependences on the ac amplitude 
and frequency at fixed temperature for arbitrary dc bi- 
ases allowing one to adjust the asymmetry of the PPP. 

The aim of this paper is to physically analyze the 
tilted-ratchet problem on the basis of the single-vortex 
model in order to determine those "intrinsic" tilted- 
ratchet effects in the vortex dynamics which arise from 
the tilt of initially symmetric WPP as the only reason. 
In addition, this model allows one to study theoretically 
the exact ratchet behavior of absorbed power at strong 
ac amplitude and arbitrary frequency, i.e., the subject 



which has not been studied in any known for us previous 
theoretical work even for usual vortex ac response. As 
a result, two groups of new findings have been obtained. 
Exact formulas for (i) the dc voltage ratchet response 
and (ii) absorbed power in ac response will be discussed 
as functions of ac current amplitude and frequency as 
well as dc current induced tilt, in a wide range of corre- 
sponding dimensionless parameters. Experimentally, the 
obtained results can be verified on superconducting films 
with a WPP, similar to those used in Refsi 14 i 15 . From the 
viewpoint of basic research, it will be pointed out which 
new ratchet effects in the vortex dynamics appear even 
within the single-vortex approximation. Besides, a fur- 
ther development of the theory towards the consideration 
of an asymmetric-potential ratchet will allow one to dis- 
tinguish tilted-ratchet effects from potential asymmetry- 
induced effects when considering an asymmetric WPP in 
the presence of a dc tilting bias. 

The organization of the paper is as follows. In Sec. II 
we introduce the model and summarize the expressions 
for dc and ac ratchet responses, obtained in terms of a 
matrix continued fraction. In Sec. Ill we graphically ana- 
lyze these quantities as functions of their driving param- 
eters, namely dc bias, ac amplitude, and frequency. In 
Sec. IV we discuss in detail two limiting cases, the adia- 
batic and nonadiabatic regimes, and explain peculiarities 
in the ratchet responses on the basis of either the static 
current- voltage characteristics or solution in terms of the 
Bessel functions, respectively. In Sec. V we conclude with 
a general discussion of our results outlining the difference 
between the intrinsic and tilted-ratchet models elucidat- 
ing their applicability and drawing parallels between our 
theoretical results and a recent experiment^. 



II. FORMULATION OF THE PROBLEM 

Our theoretical treatment of the system, schematically 
shown in Fig. 1, relies upon the Langevin equation for a 
vortex moving with velocity v in a magnetic field B = nB 
(B = |B|, n = nz, z is the unit vector in the z direction 
and n = ±1) which, neglecting the Hall effect, has the 
form 

r?v = F L + F p + F th , (1) 

where = n($o/c)j x z is the Lorentz force ($o is 
the magnetic flux quantum, and c is the speed of light), 
j = jW = j + j ac coswt, where j dc and ] ac are the 
dc and ac current density amplitudes and uj is the an- 
gular frequency. F p = — \7U p (x) is the anisotropic pin- 
ning force, U p {x) — (U p /2)(1 — coskx) is the periodic 
washboard pinning potential with k — 2ir/ a 20 23 , where 
Up is its depth and a is the period (see Fig. 2). F t h 
is the thermal fluctuation force and 77 is the vortex vis- 
cosity. We assume that the fluctuational force F t h(t) is 
represented by a Gaussian white noise, whose stochas- 
tic properties are given by the relations (F t h,i(t)) = 0, 
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FIG. 1: The system of coordinates xy with the unit vectors 
x and y is associated with the WPP channels which are par- 
allel to the vector y. The coordinate system x'y is associ- 
ated with the direction of the transport current density vec- 
tor j — j dc + j ac cos cut, a is the angle between j and y, /3 is 
the angle between the average velocity vector (v) and j. (F p ) 
is the average pinning force provided by the WPP, Fl is the 
Lorenz force for a vortex, and B is the magnetic field vector. 
Inset: a schematic sample configuration in the general case. 
Experimentally deducable values are the voltages E x and E y . 



{F th ,i(t)Fth,j(t')) = 2Tr]5ijS(t - t'), where T is the tem- 
perature in energy units, (...) means the statistical aver- 
age, Fth,i{t) with i — x or i — y is the i component of 
F t h(t), and Sij is Kronecker's delta. 

The Langevin equation (JlJ with the Hall term has 
been solved in Ref. 16 in terms of a matrix continued frac- 
tion. Neglecting the Hall effect, which is usually small in 
conventional type-II superconductors, e.g., in Nb films, 
below we summarize only the final expressions needed 
for the subsequent analysis in Sees. Ill and IV, where 
the main quantities of physical interest are (i) the time- 
independent (but frequency-dependent) dc electrical field 
response and (ii) the stationary ac response on the fre- 
quency u, independent on the initial conditions. Both 
these are determined by the appropriate components of 
the average electric field induced by the moving vortex 
system, (E(t)), whose time-independent dc components, 
(Et% and (E dc >ffi&, are 

f (E^)% = npfUf* - (sinx)£) = p fV $jf 

{ (2) 

I (£x C >0 = PfJz- 

where pf = B^o/nc 2 is the flux- flow resistivity, j c = 
cU p k/2$ , 3 d y = fcosa, j dc = fsina, f = |j dc | and 
Vq is the (u!,j dc , j ac , T)-dependent effective mobility of 
the vortex under the influence of the dimensionless gen- 
eralized moving force j dc = nj dc /j c in the x direction 



FIG. 2: Modification of the effective pinning potential U (x) = 
U p (x) — Fx with gradual increase of the Lorentz force compo- 
nent in the z-direction F, where U p (x) = (U p /2)(1 — coskx) 
is the WPP with its depth U p and period a — 2n/k. As the 
initial WPP is symmetric, i.e., U p (— x) — U p (x), it can es- 
tablish ratchet properties only in the presence of an external 
dc bias F invoking its tilt. Depending on the bias value, in 
the absence of an ac current and assuming T — K for sim- 
plicity, two qualitatively different modes in the vortex motion 
appear, (i) If F < F p , thought the initial potential well is 
tilted, it maintains the average vortex position, i.e., the vor- 
tex is in the localized state. At the critical tilt value, i.e., 
when F = F p , the right-side potential barrier disappears, (ii) 
At last, when F > F p , the vortex motion direction coincides 
with the direction of the moving force F, i.e., the vortex is in 
the running state with an oscillating instantaneous velocity 
with a frequency w U (£ d ) 2 — I (see Sec. Ill for details). 

being 

= 1 - (sinx)^. (3) 

The term (sinx)g represents the time-independent static 
average pinning force, given by Eq. (24) of Ref J^. 

The nonlinear power absorption in the ac response per 
unit volume and averaged over the period of an ac cycle 
is given in accordance with Eq. (85) of Ref J£ by the 
following expression 

P{lj) = pf/2- (j ac ) 2 [sin 2 a + cos 2 aReZi(w)]. (4) 

where 

Zi(w) = 1 - (sinx) tl /j ac (5) 

is the nonlinear impedance with the term (sinx)ti being 
the time-dependent dynamic average pinning force, de- 
termined by Eq. (24) in RefJ£, j ac = nj% c /j c , jy C = 
j a cos a, = j a sin a, and j a — |j ac |. 

III. MAIN RATCHET RESULTS 

The main goal of this section is to present results 
of a detailed theoretical study of the ratchet properties 
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FIG. 3: The ratchet voltage E d versus £ a for a set of biases £ d = 0.01; 0.1; 0.3; 0.5; 0.8; 1.05; 1.2, as indicated, in the adiabatic (a), 
intermediate (b) and high-frequency (c) regime. 



of the (dc+ac)-driven nonlinear time- and temperature- 
dependent vortex dynamics within the frames of the ex- 
act solution of the problem, presented recently in our 
work 1 -. Here we study two main physical quantities 
which can be measured experimentally for our model, 
the dc electric field E d and the ac power absorption P (in 
units of pf), as functions of their dimensionless external 
driving parameters, namely dc bias £ d = j d /j c , ampli- 
tude £ a = j a /j c , and frequency of the ac input fi = ujt 
with t = 2rj/U p k 2 being the relaxation time^. Whereas 
Eqs. ©-([5]) are written for any arbitrary angles a, in 
order to get a more simple and clear physical interpreta- 
tion of the obtained results, below we put emphasis on 
the case when a — 0°, i.e., when both the currents flow 
along the WPP channels provoking the vortex movement 
perpendicular to them. As a result, below and through- 
out the paper we consider only the y components for both 
ratchet, dc and ac responses, omitting the index y and 
(...) to simplify the notation. 

The single vortex approximation used in Ref.— sup- 
poses the WPP period, a, to be large in comparison with 
the effective magnetic field penetration depth, A, and the 
temperature low enough to prevent smearing of singular- 
ities in the ratchet responses. To accomplish this, until 
not stated otherwise, all the figure data are calculated 
for the dimensionless inverse temperature g = U P /2T = 
100i£ representing a reasonable value, experimentally 
achievable, e.g., for thin Nb films either grown on facetted 
sapphire substrates^, or furnished with nano-fabricated 
PPP landscapes", where U p ~ 1000 4- 5000 K and 
T«8K. Assuming a triangular vortex lattice matching 
the PPP landscape at a magnetic field B w 10 mT, the 
pinning structure's period is a ss 400 nmi£. For a Nb 
film with a thickness d ~ 60 nm, A(0) ~ 100 nmi£ (de- 
pends on temperature and the film's quality), so that the 
condition d < X < a can be experimentally satisfied. 

Two groups of our findings commented below refer to 
the amplitude and frequency dependencies of the ratchet 



responses E d and P. The dependencies are complemen- 
tary to each other allowing one to describe quantitatively 
E d and P in the whole (£ a , f2, £ d )-space, as detailed next. 



A. Electric field dc response 

One of the main questions in the study of the function 
E d (t; a \f; d ,n) given by Eq. (J2)) relies upon the determi- 
nation of the frequency and dc bias dependences of the 
ac amplitude threshold value, £?(f2,£ d ), which can be 
considered as an ac critical current magnitude for the dc 
ratchet response E d such that E d = for £ a < . To ac- 
complish this, we begin the graphical analysis with the ac 
amplitude dependence of the dc ratchet response consid- 
ering specific features in E d (£ a \£l, £ d ) at low (ft — 0.01), 
intermediate (f2 = 1), and high (ft = 3) frequencies for 
similar tilts, as depicted in Fig. 3. 

Consider at first the curves in Fig. 3. a with SI = 0.01 
which corresponds to the adiabatic ratchet response. At 
large £ a , for all the curves tt < 1) ~ £ d . These 

values will be explained in Sec. IV in a rather simple man- 
ner as this asymptotic behavior follows from Eq. (jl"5j) for 
E + (£, d + £ a coswi) at £ a -> oo with (+) denoting the 
even component of E regarding the change £ Q — > — £ a . 
At small £ a values, we observe different behavior for 
curves with £ d > 1 and £ d < 1, namely for £f = 1.05 
and £ d = 1.2 the ratchet response is a threshold- free 
one, whereas a threshold value, separates the non- 
dissipative and dissipative states at £ d < 1. The magni- 
tude of the threshold is a decreasing function of £ d and, 
in fact, is equal to = 1 — £ d which is evident for the 
adiabatic case. The physical reason of the above differ- 
ence follows from the fact that at £ a = and £ d > 1 the 
vortex is in the running state with a slightly oscillating 
instantaneous velocity dx/dt and thus, nonzero electric 
field E, whereas for £ d < 1 the vortex is localized in one 
of the WPP wells. 
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FIG. 4: The frequency dependence of for the dc- 
tilted cosine pinning potential at different biases £ d = 
0.1; 0.3; 0.5; 0.65; 0.8;, as indicated. The navy and orange 
(online) dash lines represent rough separations between low 
(fi <C 1), intermediate (O ~ 1), and high frequency (O ^> 1) 
regimes. Inset: the nonlinear transition from the adiabatic to 
high-frequency regime in detail. The curves behave qualita- 
tively similar to those obtained experimentally in Ref.— on 
superconducting Pb films with a non-tilted ratchet pinning 
potential. 



Transferring from low (f2 = 0.01) to intermediate 
(51 = 1) frequencies, several new distinctive features ap- 
pear in Fig. 3.b in comparison with the adiabatic case. 
First, the threshold values = l,£ d ) are larger than 
those ensuing for SI = 0.01 at similar subcritical dc tilts, 
i.e., £ d < 1. To illustrate this in detail, we plot the £™ (fl) 
dependence in Fig. 4 for a set of biases. All the curves 
demonstrate qualitatively similar behavior, i.e., a zero 
plateau at £ a < a linear dependence at large £ a > , 
and a nonlinear transition in between at £ a > £™ . These 
segments correspond to the adiabatic, intermediate, and 
high-frequency modes which are roughly separated by the 
straight lines O ~ 0.1 and ft ~ 1, respectively. It should 
be noted, that the curves £"(51) in Fig. 4, calculated in 
the present work for the dc-tilted cosine pinning poten- 
tial, are qualitatively similar to those obtained experi- 
mentally on superconducting Pb films with a non-tilted 
ratchet PPP (see RefJ^ and Fig. 5 therein). The tran- 
sition frequency from the adiabatic to nonadiabatic case 
has been found at about 1 MHz for that system-^. 

Second, a difference in the E (£ a |£ , O) behavior ap- 
pears between 0.4 < ^ idd i e ~ 0.7, which looks like 
damped oscillating curves, and the curves at £ d < 0.4 
and £ d > 0.7, which look like curves with phase-locked 
regions (steps) in £ a . Whereas at small £ d phase-locked 
regions ensue at E d = 0, at strong biases £ d > 0.7 
these flat segments appear at E d = 1. In Sec. IV this 
will be discussed in detail within an approximate Bessel- 
function approach, originally presented in Refs i 25 i 26 and 
used later— for O = 1. In fact, the authors of Ref. 27 



numerically calculated the Langevin equation in order to 
obtain the so-called dynamical current-voltage character- 
istics (CVCs) of the resistively shunted Josephson junc- 
tion model, which is equivalent to the model used here 
to analytically derive the dc ratchet-response solution 
E d (^ a \£_ d ,il) (notice the curve E d (£ a \£ d = 0.8, Q = 1) 
in Fig. 7. a of Refill). In our model, the corresponding 
Langevin equation with noise has been exactly solved in 
RefJ£ (see Eq. (9) therein) and will be discussed in the 
noise-less limit (Eq. (|22p in Sec. Ill) to clarify the origin 
of the steps in Fig. 3.b. 

Lastly, consider the curves in Fig. 3.c at O = 3. From 
the above discussion it is clear that for £ d > 1, the ratchet 
responses E d are continuously oscillating curves without 
thresholds. On the contrary, for £ d < 1 the responses 
have thresholds whose magnitudes decrease with increas- 
ing £ d . As in the previous case, an interesting property 
of the dependence E d (£ a \!; d , ft) in Figs. 2.b,c is the possi- 
bility for E d to decrease periodically (sometimes down to 
zero) with increase of the driving amplitude £°. Such a 
behavior of E d is in contrast to the behavior of the usual 
dc-driven CVCs, even though in the presence of £ a , as 
these are always increasing functions of £ d . 

Now we turn to the analysis of the frequency depen- 
dences of the ratchet response E d taken at fixed £ a and 
£ d , as represented in Fig. 5 for small (£ a = 0.5), inter- 
mediate (£ a = 1) and strong (£ a = 3) ac amplitudes. 

At small ac drives £ a = 0.5 (see Fig. 5. a) the curves 
vanish regardless of the frequency at small tilts £ d < 0.5. 
This behavior is evident, since if both the tilt value and 
ac drive amplitude are small, the vortices are localized 
at the bottoms of the WPP wells which results in a non- 
dissipative state. With the gradual increase of the bias, 
for determinacy from £ d = 0.5 to £ d = 0.8, the situa- 
tion changes dramatically. At low frequencies the volt- 
age drop gets substantially higher, whereas a zero- voltage 
tail spreads over the high-frequency range. The former is 
a consequence of the running vortex state, whereas the 
latter is a clear signature of the localized vortex state. 
These regions are separated by a threshold frequency 
il c = Sl(£ d , which is in fact already depicted in Fig. 4 
and, as evident from the latter plot, is strongly dependent 
on both, £ d and £ a (compare with the curves in Fig. 5.b 
for £ a = 1). The tilt, £ d , determines the asymmetry of 
the WPP and the time needed for a vortex to get from 
one to the next WPP well, whereas £ a represents the 
ac driving force for a vortex which also competes with 
the height of the initially symmetric WPP. This physi- 
cally means that, if the ac driving frequency is much 
less than the depinning frequency ui p ~ t -1 , the running 
state of the vortex appears and it can visit several poten- 
tial wells during the ac period. For a fixed ac amplitude 
£ a and frequency O, the number of visited wells increases 
strongly with the increase of the tilt, thus resulting in a 
shift of the threshold frequency towards higher tt. An- 
other interesting feature in the E d (il) curves appears as 
a maximum at il ~ 1. Its magnitude increases with in- 
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FIG. 5: The voltage E d versus 51 for a set of biases £ d — 0.01; 0.1; 0.3; 0.5; 0.8; 1.05; 1.2, as indicated, at small (a), intermediate 
(b) and strong (c) ac drives. 



crease of the frequency. 

The behavior of f) plotted for £ a = 3 in 

Fig. 5.c, as representative for strong ac drives, can be 
summarized as follows. In the adiabatic limit, when 
51 < 1, the function E d (Q\^ d , £ a ) coincides with the tilt 
values £ d . At high frequencies 51 > 1, the curves either 
attenuate rapidly for subcritical tilts £ d < 1 or approach 
a constant value for £ rf > 1. In the intermediate regime, 
when 51 ~ 1, the curves oscillate until the maximum is 
reached, followed by the subsequent rapid decrease of the 
function i? d (51|£ d , £ a ). Physically, this corresponds to the 
pronounced reduction of the time T u /2 with T u being 
the ac period, over which the moving force on the vor- 
tex keeps its direction. As a result, the vortex may not 
longer visit several WPP wells, since the vortex displace- 
ment during T w is smaller than the WPP period a, even 
at strong ac amplitudes £ a . 

Summarizing, the calculated ratchet behavior of 
E d (£ a \£ d , 51) differs substantially for a wide range of fre- 
quencies 51, ac amplitudes £ a and dc biases £ d . Through- 
out the frequency range this include a cut-off filter be- 
havior of the function E d (£ > a \t; d 1 51) with decreasing £ a 
and increasing 51. In addition, at £ d < 0.4 and £ d > 0.7 
the curves demonstrate phase-locked peculiarities remi- 
niscent of Shapiro stepsi^ as well as a damped oscillatory 
behavior at 0.4 < £ d 



t d < 7 

^middle ~ w * 1 * 



B. Power absorption in ac response 

In this subsection we first consider the behavior of the 
absorbed ac power V, calculated as function of the di- 
mensionless dc density £ d , ac density £ a , and frequency 
51 in Sec. V.E of RefJ^. Then we consider the graphs 
for E d (£ a \Z d ,fl) and ReZi(£°|£ d , 51) in comparison with 
each other. 

Before entering the discussion it is useful to remem- 
ber^ that 51 = lot, where f is the relaxation time and 
ujp ~ f" 1 is the depinning frequency calculated in the 



linear ac approximation at £ d = 0. We have to point out 
that the bias dependences for the critical frequency 51 c 
introduced in the previous subsection and the depinning 
frequency to p are qualitatively opposite. Whereas the 
frequency dependence of the first Shapiro-like anomaly 
in E d (^ d ) is given by the well-known expression 51 = 
V (£ d ) 2 — l 24 f° r overcritical biases, the dependence for 
the depinning frequency w p (£ d ) = u p yl — (£ d ) 2 has been 
obtained in the linear ac approximation at £ a — > and 
for subcritical tilts at T = only recently 28 -. Thereby, in 
Refi2£ it was shown that decreases with increas- 

ing £ d at < £ d < 1. The physical meaning of w p (£ d ) 
at low temperature (g 3> 1), i.e., when the creep factor 
voo (see Eq. (107) ir>i&) is small, follows from the fact 
that for low frequencies ui <C w p (£ d ) (or 51 <C 1) pinning 
forces dominate and the vortex resistivity response p v , 
being proportional to the absorbed power V, is nondissi- 
pative (see Eq. (2) in Ref J£), whereas at high frequencies 
u ^> w p (£ d ) (or 51 3> 1) frictional forces dominate and p v 
is dissipative. The power absorbed per unit volume and 
averaged over the period of an ac cycle V{lo) was calcu- 
lated ir>i£ (see Eq. (84) therein), and can be written as 



P(u } ) = (p f /2)[(£) 2 + (a 2 ReZ 1 (u)} ) 



(6) 



where Zi(£ a , £ d , 51, a, g) is the nonlinear frequency- and 
dc and ac amplitude-dependent impedance. In RefJ£ it 
was shown that Z\ plays the same role for the ac response 
as Vq for the dc response. 

Proceeding now to the analysis of the dependences 
V{£ a \£ d , 51), let us recall that to accomplish this in the 
case a — 0°, it is sufficient to calculate the ac resistiv- 
ity pi = ReZi(£ a |£ d , 51) as function of its parameters, as 
presented in Figs. 5 and 6. 

As for the dependences E d (^ a \^ d , 51) depicted in Fig. 4 
and corresponding to low, intermediate, and high fre- 
quencies, the curves pi(£ a |£ d , 51) are plotted at the same 
£ d values. Consider at first the curves Pi(t; a \t; d , 51) in 
Fig. 6. a at 51 = 0.01 which correspond to the adiabatic 
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FIG. 6: p\ versus £ a for a set of frequencies Q and biases £ d = 0.01; 0.1; 0.3; 0.5; 0.8; 1.05; 1.2, as indicated, in adiabatic (a), 
intermediate (b) and high-frequency (c) regime. 



case w <C w p (£ a ). Here localized vortex states ensue 
at £ d < 1, and running states appear at £ d > 1. For 
small tilts, £ d < 1, an absorption threshold appears in 
the pi(£ a |£ rf ,ft = 0.01) curves at £ a = Q. Here £? coin- 
cides with the critical ac magnitude for the correspond- 
ing curves E d (£, a \£, d < l,ft = 0.01). The physics of this 
threshold was earlier discussed for ft = 0.1 (see RefJ£, 
Fig. 6) and for our adiabatic case (ft = 0.01) may be 
connected with the u! p (£ d ) dependence. At subcritical 
ac drives and overcritical biases, e.g., at £ d = 1.05, as 
a consequence of the running vortex state, pi acquires 
large values. With the gradual increase of £ a the curves 
exhibit a weak minimum and finally approach unity. 

In the case of intermediate frequency (see Fig. 6.b for 
Q = 1) the curves /di(£ a |£ d , ft) start from a nonzero value 
regardless the bias £ d . A different behavior at subcriti- 
cal and overcritical tilts should be noted. At overcritical 
biases the running state for the vortex appears and the 
response is a consequence of this motion with a slightly 
oscillating instantaneous velocity. On the contrary, at 
subcritical tilts the vortex is in the oscillating state and 
if the frequency of the external excitation exceeds the 
depinning frequency w p (^ d ), a nonzero response appears 
as a result of the averaging of this oscillations. The bias 
dependence at £ d < 1 in the limit of small ac drives and 
ft <C 1 has been discussed above. Proceeding with the 
analysis of the curves shown in Fig. 6.b it should be noted 
that at small biases £ d < 0.5 the curves increase mono- 
tonically, whereas at sufficiently large tilt values, e.g., 
at £ d — 0.8; 1.05; 1.2 a pronounced minimum appears. 
At overcritical biases this can lead to a sign change in 
/9i(£ a |£ d ,ft) in both £ a and (, d . This evidence is repre- 
sented in detail in Fig. 7 overcritical tilts and weak ac 
drives. Figure 7 illustrates the presence of the singular 
point £ d = Vl + ft 2 in pi(£ d ,t; a ) at nonzero both the ac 
current amplitude and the temperature. A similar singu- 
larity in the Josephson junction impedance problem has 
been discussed previously in the limit of small ac am- 



plitude at zero temperature 2 ^. Here we show that such 
a singularity does not vanish at nonzero temperature in 
the presence of subcritical ac amplitude £ a > 0. This evi- 
dence was left out of the scope of that work 2 ^. Proceeding 
with the description of this anomaly it should be noted, 
that the dependence of pi(£, a \£, d on £ a is less sharp 
than in £ d . As far as high temperatures smear the singu- 
larity in both £ a and £ d , the minimum can be more pro- 
nounced if considered at lower temperatures (g > 100). 

Considering /Oi(£ a |£ d , ft) at high frequencies in Fig. 6.c, 
it is evident that /9i(£ a |£ d ,ft) approaches unity even at 
small £ a which again corresponds to the case w > uj p (£, d ). 




P1 



FIG. 7: The ac resistivity p± versus C and £ d at Q = 1 
demonstrating a sharp singularity in £ d . This results in a sign 
change in the function pi(£ d ,Q) at overcritical tilts £ d > 1. 
The position of the minimum can be quantitatively calculated 
as £ d = y/l + Q. 2 . 
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FIG. 8: pi versus for a set of biases £ d = 0.01; 0.1; 0.3; 0.5; 0.8; 1.05; 1.2, as indicated, at small (a), intermediate (b) and 
strong (c) ac drives. 



We go on with the analysis of the frequency depen- 
dence of pi, represented in Fig. 8 for small £ a = 0.5, 
intermediate £° ~ 1 and strong £ a = 3 ac drives for the 
same set of biases £ d . 

At small ac drives (see Fig. 8. a) the curves demonstrate 
either monotonic behavior for £ d < 0.5 or pronounced 
non-monotonic behavior for £ d > 0.5. The monotonic 
curves at £ d — agree with the results of Coffey and 
ClemSi who calculated in linear approximation in £ a the 
temperature dependence of the dcpinning frequency in 
a nontilted cosine pinning potential. In contrast to this 
monotonic behavior, the nonmonotonic curves (£ d > 0.5) 
demonstrate two characteristic features. First, a pro- 
nounced power absorption in the adiabatic regime. Sec- 
ond, a deep minimum in the power absorption at fi ~ 1. 
The appearance of this frequency- and temperature de- 
pendent minimum was discussed in more detail in Ref 
With the gradual increase of £ a , the value p\ at low fre- 
quencies remains the same, whereas the minimum shifts 
towards higher frequencies (see Fig. 8.b). In addition, at 
frequencies 17 ~ 0.5 peculiarities in the curves become 
more pronounced. These can be smeared in turn when 
considered at higher temperatures (g « 100). 

At strong ac drives, as represented in Fig. 8.c for 
£ a = 3, already at very low frequencies ft <C 1 all the 
curves p\ acquire large values and approach unity at 
high frequencies Q 3> 1. Even though peculiarities in 
the dependence pi(fi) seem to be pronounced in the in- 
termediate frequency range, they are sufficiently weak 
in comparison with those in Fig. 8.b when £ a = 1. A 
further increase of £ a leads to slightly distorted curves 
throughout the frequency range. 

The main results of this subsection can be summarized 
as follows. The power absorption in the ac response has 
been considered in terms of the ac resistivity pi as a func- 
tion of its driving parameters, £ a , ft, and £ d . While in the 
limiting cases of small ac current in the absence of dc bias 
the well-known results of Coffey and Clem follow^ and 
at strong dc biases a large power absorption results, in 



agreement with the curves reported previously^, the ap- 
pearance of a sign change in pi(£ d , ft) at a certain range 
of ac drives £ a at overcritical biases t; d is predicted for 
the first time in the present work. 



IV. DISCUSSION 

In this section, two physically different limiting cases 
at low and high frequencies will be considered at zero 
temperature (i.e. g — ¥ oo) to augment the previous anal- 
ysis with a more intuitive and visual interpretation. 

The first case we consider is the adiabatic regime with 
SI <C 1. To discuss the dc ratchet response E in this limit, 
we employ static CVCs (see Eqs. © and (ITU)) below) with 
j = £ a cos uit, jo = £ d and average E + (j + jo) over the 
driving period T u = 2it/uj. In this limit (w — > 0), while 
£ a > 1 the vortex may visit many potential wells of the 
WPP during the time i.e., when the moving force 

on the vortex keeps its direction. 

The second case is the non-adiabatic case with S7 3> 1 . 
In this limit, the vortex displacement during the time T u 
may be smaller than the WPP period a, even at strong 
ac densities £ a and the Langevin equation for a vortex 
can be solved in terms of the Bessel functions as detailed 
below. 



1. Adiabatic case 

We consider the vortex motion in a WPP, tilted by a 
dc bias £ d < 1, and subjected to an ac drive with u — s- 0. 
Our objective is to qualitatively point out why a rectified 
dc voltage appears in response to the ac input £ a . To 
accomplish this, we first consider the dc CVCs for the 
cosine WPP (see Eqs. ©) in the absence of an initial 
tilt, viz., 

E (j) = jv (j), (7) 
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FIG. 9: The functions vo(j) (a) and Eo(j) (b) in the absence of a tilt in the adiabatic case at zero temperature, (c) In the 
presence of a dc bias jo = 0.8, E(j + jo) (solid) and E(—j + jo) (dash) result in the appearance of even E + (j + jo) (dash dot) 
and odd E~ (j + jo) (dot) components (see Eq. 1)120 ). The horizontal and diagonal thin dot lines represent asymptotics E = 0.8 
and E — j, respectively, (d) E + (j + jo) (dash dot) is plotted together with the ratchet response E r (crosses, Eq. (|19p ) and the 
exact solution for E d (j,Q = 0.001, jo = 0.8, g = 1000) (solid, Eq. ©). 



where — oo < j < +oo and 



and 



^o(j') 



0. 



lil > i, 
lil < l, 



(8) 



is the well-known nonlinear vortex mobility under the in- 
fluence of the dimensionless generalized moving force j in 
the ^-direction (see Eq. 52 in Ref. 1 ^). From Eq. (JSJ fol- 
lows, that fo(J) is an even function of j, i. e. vq{J) — 
Vo{~j), whereas Eo(J) is an odd function of j, i. e. 
E o{j) = -E (-j). The functions v (j) and E (j) are 
shown in Fig. 9. a and Fig. 9.b, respectively. 

As the cosine WPP is symmetric, i. e. U p (x) — 
Up{~x) 1 it can establish ratchet properties only when 
tilted, i.e., we change j — > j + jo, where jo is the tilting 
dc bias. It is easy to see that depending on the sign of 
jo, the tilt leads to the shift of v§ and E along the j-axis 
by the value |jo| to the left (for j > 0) or to the right 
(for jo < 0). This is illustrated in Fig. 9.c where jo = 0.8 
for definiteness. In this case 



where 



Eoti) -> E(j+j ) = {j+jo)v(j+jo), 



Vi-VO' + jo) 2 , 



(9) 



j>l-3o, 

Mj+h) = < 0, -1 - j < j < 1 - j , 

I Vl-l/O'+J'o) 2 , 3 < (-1 - Jo), 

(10) 

Whereas Eo(j) and vo{j) are odd and even functions of 
j, respectively (see Fig. 9), from Eqs. (|9]) and (ITUl) follows 
that in the presence of a tilt E(J + jo) and v(j + jo) are 
neither even, nor odd in j for jo 7^ 0. In the following it 
is suitable to present E(j ± jo) as 



E(j ± j ) = E+(j ± j ) + E-(j ± j ), 



(11) 



where 



E ± (j ± jo) = [E(j ± j ) ± E(-j ± j )]/2 (12) 



E ± (j±j )=±E ± (-jT3o) 



(13) 



are the even and odd parts of E(j ± jo) with respect to 
change j -> - j. 

On the other hand, it is clear that 



E(j±j ) = -E(jT3o), 



(14) 



because when we change the sign of the full current [i. e. 
(j ± jo) — » (—j T jo)], then E will be an odd function 
of the full current. If we then apply the i-operation to 
E given by Eq. (|14[) and take into account Eq. (fT2|) . we 
arrive at the important conclusion that 



E ± (j ± JO )= T E ± (jT Jo)- 



(15) 



From Eq. (JTSJ) it follows that E + (j+j ) = -E + (j~j ). 
This means E + (j + jo) changes its sign when changing 
the sign of jo, whereas E~(j — jo) = E~(j + jo) does not 
change its sign. From the physical viewpoint it means 
that E + (j + jo), even in j and odd in jo, is responsible 
for the ratchet response, whereas E~ (j+jo), which is odd 
in j and even in j , describes the usual CVCs response, 
analogous at j 3> 1 to that at jo = 0. Actually, if we take 
into account that from Eq. (fTOjl follows lim^oo v{±j + 
jo) = 1, then from Eq. (fT3|) at once follows that 



.lim E+(j+j ) 

j->±oo 



lim E (j 

— >±oo 



-Jo) 



Jo, 



J- 



(16) 



(17) 



At last, let us perform the change j — > j cos cot and con- 
sider the function E(j cos ut + jo), since this represents 
a more close correlation with the exact results obtained 
by using Eq. @. To derive the average dc ratchet solu- 
tion E r in response to to the input current j cos tot + jo, 



10 



one needs to integrate the function E over the ac current 
period 



€2 t% 



E r = — dt E(j cos u>t + j Q ). 

J-u) Jo 



(18) 



Equation ([15} can be reduced to the sum of two integrals 



7T 



tt/2 pir/2 

dtp E(j cos tp+j )+ / dip E(-j cos ip+j )}, 
Jo 

with the integrals to be taken over tp — ujt only where 
p(j cos ujt + jo) is nonzero in accordance with Eq. (1101) . 
i.e., 



JO - 



<p > 



arccos( — 
arccos( ~ T J ° ' 



J > i- jo, 
j < -1 - jo, 



(20) 



Equation (1191) can be rewritten in another equivalent 
form 



E r 



tt/2 



d(/3 E+(jcosip + j Q ), 



(21) 



which represents the even component of the dc ratchet re- 
sponse and will be compared next with the temperature- 
dependent response E d (^ a ). 

We now take a closer look at Figure 9 considering 
firstly the curves E(j +j ) (solid) and E(—j+j ) (dash) 
in Fig. 9.c. As it follows from Eq. (fl6l) . E(j + jo) is zero 
at — 1 — jo < j < 1 — jo as far as E(—j + jo) vanishes 
at — 1 + jo < j < 1 + jo . A rapid increase in both the 
functions in j in the vicinity to ±(1 — jo) an d ±(1 + jo) 
should be noted. In the adiabatic limit (see Eq. (|12l) ). 
E + (j +jo) and E~(j + jo) are shown in Fig. 9.c by dash 
dot and dot lines, respectively. They inherit both the 
bump-like peculiarities in E(±j + jo) at ±(1 + jo)- Next 
we turn to Fig. 9.d where adiabatic (Eq. (fl6|) . dash dot), 
approximate (Eq. (fT9|) . crosses), and exact (Eq. ([2]), solid) 
ratchet responses are shown together. The asymptotic 
behavior of all the curves is in agreement with Eq. . 
Peculiarities are highly pronounced in the adiabatic so- 
lution (fT2j) . calculated in fact in response to the square 
wave ac current. In contrast to this, the approximate 
ratchet solution (fT9|) for the cosine ac current practically 
coincides with the exact solution ([2]) calculated in the 
limit of very small frequencies ui = 0.001 and very low 
temperature g = 1000. 

By this way, we conclude that the simple approach in 
terms of the static CVCs can explain qualitatively the 
form of the curves E d (C) in Fig. 3 .a in the adiabatic 
limit. A good quantitative agreement between the ap- 
proximate and exact solutions is revealed at very low 
frequencies and temperatures. 



2. Nonadiabatic case 

In order to explain qualitatively the results of the 
exact calculations of E d (£ a \£ d , f2) at low temperatures 
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FIG. 10: The functions J (C) and Ji(? a ) (dash) are plotted 
together with straight lines £ d = 1.2 and £ d - fi = ±0.2. 
The phase locked regions in E d (£ a \M = l,C d = 1-2,5 = 100) 
(blue line) are connected by vertical lines for clarity. The 
light segments of the stripe on top of the plot indicate phase 
locked regions with £° denoting the roots of Eq. (I27[) for n=0 
and n=l. A better fit of the noiseless approximate solution 
(23) can be achieved if the ratchet response ([2]) is calculated 
at very low temperatures, such as for g = 1000 (black line). 



(g = 100), presented in Fig. 3b, c for intermediate and 
high frequencies (f2 ~ 1 and f2 3> 1) of the ac driving 
force, we use a more simple approach which ignores the 
noise (g —> 00). In this limit, the equation of motion for 
the dimensionless vortex coordinate x (see Eq. (9) in^) 
reduces to 



dx/dt + sinx = £ d + £ a cos Vtt 



(22) 



which is analogous to the well-known equation of mo- 
tion for the phase difference in the ac-driven resistively 
shunted Josephson-junction model^i at zero tempera- 
ture. For high and intermediate frequencies (f2 3> 1 and 
fl ru 1) and for £ d > 0.5 Eq. (|2"2"|) can be approximately 
analyzed and solved in the spirit of the ansatz of Refi^. 
In these limiting cases we simply assume that the velocity 
of the vortex is sinusoidal in accordance with 



dx/dt = (dx/dt) + C cos nt 



(23) 



and determine the constant (dx/dt) by requiring that 
dx(t)/dt satisfies Eq. ([22]). Having integrated Eq. (j23|) 
we obtain 



i(t) =x + (dx/dt)t + (£°/fi) sinOt, 



(24) 



where xo is a second constant to be determined. Substi- 
tuting Eq. (pM} into Eq. (f2"2"]) and using the expansion of 
sin x in a harmonic series according to formulas originally 
suggested by Shapiro et. al. in Ref,— , i.e., 

cos(xsin<^) = J2kL-oo Jk(x-)cos(ktp), 



sin(xsin</?) = J2kL-oo ^fc(x) sin(fc</?) 



(25) 
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FIG. 11: The functions Jo(0 and Ji{C) (dash) are plotted 
together with straight lines ±£ d (= 0.1) and Q - £ d = 0.9 to 
illustrate the phase locked regions in E d (£ a \Sl = 1) at £ d = 
0.1. Other details are similar to those in Fig. 10. 



where Jk (x) is the fc-th order Bessel function, we obtain 

oo 

(dx/dt) = £ d - J k(Z d /ty sin[x + ((dx/dt) + fcfi)t] 

k— — oo 

(26) 

from which (dx/dt) can be found self-consistently. In 
Eq. (|26l) xo is an arbitrary coordinate. For values (dx/dt) 
that are not integral multiples of f2, the term sinx does 
not contribute a constant component to Eq. (|22p. and 
equating the constant terms yields (after averaging over 
one period for locking into the n-th region) 



S, d - rift = J„(f /ti)(-l) n sinx . 



(27) 



In order to use these results for the explanation of the 
behavior of E d (£ i a \£ i d = 1) (see Fig. 3.b) we consider 
two cases. The first case considers £ d = 1.2 for which at 
£ Q = the vortex is in the running state. From Fig. (9) 
we see, however, that for < £ a < the value of E d is 
locked into the phase of the ac periodic driving at n = 1 
(see Eq. (f27|l). From Eq. (gZJ) follows that for n = at 
£ d = 1.2 and S7 = 1 this equation has no solution at any 
value of xq. For n = 1 we have 



Ji(f)sinx 



-0.2 



(28) 



at values which satisfy the condition Ji(^ a ) = ±0.2 
and strictly correspond to the phase-locked regions for 
this curve. Between these regions E has a bump-like form 
with an increasing width and decreasing height with the 
increase of For n = 2 Eq. (|2T|) has no solution. 

The second case deals with £ d = 0.1, for which at 
£ a < Q the vortex is in the locked state with n = 0. In 
this case the equation Jo(£ a ) = ±0.1 has many solutions 
(see Fig. 11). Considering n = 1 leads to the equation 
Ji(£°) sinxo = ±0.9 which has no solution, however. 



V. CONCLUSION 

In this work we proposed an exactly solvable two- 
dimensional model structure for the study of the 
frequency-dependent ratchet effect in superconducting 
film with a symmetric planar pinning potential, tilted by 
a dc bias, also known as a tilted ratchet. We have the- 
oretically examined the strongly nonlinear nonadiabatic 
tilted ratchet behavior of the two-dimensional vortex sys- 
tem of a superconductor as a function of the (ac+dc) 
transport current density j, the frequency uj, and the tem- 
perature T. The nonlinear (in j) resistive behavior of the 
anisotropic vortex ensemble is be caused by the presence 
of anisotropic pinning with the symmetry of the PPP. It 
is physically obvious that such a pinning at low enough 
temperatures leads to anisotropy of the vortex dynamics 
since it is much easier for vortices to move along the pin- 
ning channels (the guiding effect in the flux-flow regime, 
which is linear in the current) than in the perpendicular 
direction, where it is necessary for them to overcome the 
pinning potential barriers. The latter is also a source of 
nonlincarity of the dc+ac responses. If under variation of 
one of the "external" driving parameters j, T, and a, the 
intensity of the manifestation of the indicated nonlinear- 
ity is weakened, this weakening will lead to an "effective 
isotropization" of the vortex dynamics, i.e., to a conver- 
gence (and in the limit of the absence of nonlinearity, to 
coincidence) of the directions of the mean velocity vector 
of the vortices and the Lorentz forced. 

It is physically clear that current, temperature, and 
angle a have qualitatively different effects on the weak- 
ening of the pinning and the corresponding transition 
from anisotropic vortex dynamics to isotropic. With the 
growth of j the Lorentz force F l grows and the height 
of the potential barrier decreases, so for j jcri,jcr2 
these barriers essentially disappear. Here j C ri,2 are the 
crossover currents for these transitions to occur regard- 
ing the right- and left-hand PPP barriers. The quantities 
jcri,2 depend on a by virtue of the fact that the prob- 
ability of overcoming the barrier is governed not by the 
magnitude of the force Fl, but only by its transverse 
component i^cosa, so that j cr i^{a) = j C ri,2(0)/cosa 
grows with increasing of a. Since an increase in tem- 
perature T always increases the probability of overcom- 
ing the pinning barrier, the transition to isotropization of 
the vortex dynamics is much steeper in T, the smaller the 
pinning barrier is. However, although general formulas 
for the ratchet responses (see Eqs. (O-©) include both, 
the angle and temperature dependences, in the present 
work we used a = 0° and g = U P /2T = 100 for simplicity. 

Proceeding now to a short description of the main the- 
oretical results, we note here that an exact analytical 
representation of the nonlinear ac-driven rachet response 
of the investigated system in terms of a matrix continued 
fraction was possible thanks to the use of a simple but 
physically realistic model of anisotropic pinning with a 
tilted cosine WPP. The exact solution obtained made it 
possible for the first time to consistently analyze not only 



12 



the qualitatively clear vortex dynamics of the adiabatic 
ratchet effect, but also the nontrivial ratchet behavior 
at intermediate and high frequencies of the ac drive. Be- 
low we turn to a short presentation of our results taken in 
comparison with experimental results, presented recently 
inRef^. 

First of all, simple inspection of our exact expres- 
sions ([5} for the dc ratchet response shows that a 
magnetic field inversion does not change the sign of 
£ d (^|f,w), as on the right side of Eqs. © the E^ y - 
components do not depend on the index n, which deter- 
mines the B-inversion. On the other hand, we should 
point out that for the adiabatic ratchet studied in detail 
in Ref.— for the asymmetric PPP without dc bias, the dc 
response changes its sign after B-inversion (see Eq. (16) 
therein). Since a clear sign change in Vdc has been exper- 
imentally observed at field inversion^ 3 -, we conclude that 
the dc ratchet response in that work should be described 
by a model with an asymmetry of the PPP—. 

An another interesting difference between the tilted 
ratchet described in the present paper and the asymmet- 
ric ratchet without a tilt"* consists in their asymptotic 
behavior at j — ¥ ±oo in the adiabatic regime. In this 
limit, the tilted ratchet adiabatic response is finite and 
equal to the tilt value, as one can see from Eq. (16) at 
jo < 1, whereas for the asymmetric ratchet this response 
is zero because lim^-too v~{j) ~ — > fA 

Comparing our frequency-dependent results with anal- 
ogous experimental findings presented in Ref^ we 
should underline that in spite of the simple WPP and 
the single-vortex approximation used in our theoretical 
model we, however, can qualitatively explain from one 
and the same point of view main experimental results 
of that work. In particular, our expressions © and dU 
describe (i) the critical ac current dependence in a wide 
frequency range covering the transition from adiabatic 
to nonadiabatic, with both, a frequency- independent 
plateau at low frequencies, a direct dependence at high 
frequencies and a nonlinear transition in between, (ii) the 
appearance of phase-locking regions in the dependence 
of Vdc on I r f, (iii) a weakening of the ratchet effect at 
extremely high frequencies, and (iv) the possibility of a 
sign change for the absorbed power in ac response within 
a certain range of the driving parameters. This is in con- 
tradistinction to different explanations used ini 3 -, which 
is a consequence of the absence of a well-defined theoreti- 
cal model due to the complexity of the PPP employed. In 



particular, due to this reason the authors^ 3 compelled to 
employ different approaches for the explanation of their 
experimental results, such as introducing sometimes the 
vortex mass, or appealing to the vortex- vortex interac- 
tion, or not taking into account the tilting parameter, or 
leaving uncommented the high frequency power absorp- 
tion behavior. 

It should be stressed that the single vortex approxima- 
tion used in this work may only be valid at small mag- 
netic fields preferably less than the first matching field, 
so no collective effects are captured in the model consid- 
ered here. Whereas vortex ratchet reversals as a function 
of field have been studied by authors in a number of sim- 
ulations and experimental works and have been explain- 
ing as a result of collective effects, such as vortex- vortex 
interactions (see, e.g., Ref. 3 — and references therein), a 
remark on the extension of our theoretical study is now 
to be made. As far as the model described in the present 
paper refers to the tilted-potential ratchet, a thorough 
theoretical description of the rocking-ratchet response in 
superconducting films with an asymmetric WPP is cur- 
rently under way and will be reported in a forthcom- 
ing publication 3 !. There will be shown that the single- 
vortex approximation can also lead to the ratchet rever- 
sals, when we consider an asymmetric WPP in the pres- 
ence of a tilting bias. Finally, we would like to stress, 
that our exactly solvable single-vortex model explicitly 
shows that many important and interesting nonlinear 
ratchet effects, which can be observed in particular at 
high frequencies, follow even from such a simple model 
for one vortex in periodic ratchet WPP. Though experi- 
mental verification of the predictions of both the models 
can be performed, for instance, on Nb thin films with 
nano-fabricated WPP landscape a 14 i 15 , the first portion 
of ratchet data still remains to be seen. 
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